Transport control in deterministic ratchet system 



Woo-Sik Son, 1,2 'EI Jung- Wan Ryu, 3 Dong-Uk Hwang, 4 Soo- Young Lee, 2 Young-Jai Park/'Q and Chil-Min Kim 2 [l 

1 Department of Physics, Sogang University, Seoul 121-742, Korea 
^National Creative Research Initiative Center for Controlling Optical Chaos, Pai-Chai University, Daejeon 302-735, Korea 
s Department of Physics, Pusan National University, Busan 609-735, Korea 
^Department of Biomedical Engineering, University of Florida, Florida 32611-6131 

(Dated: February 8, 2008) 

We study the control of transport properties in a deterministic inertia ratchet system via the 
extended delay feedback method. A chaotic current of a deterministic inertia ratchet system is 
controlled to a regular current by stabilizing unstable periodic orbits embedded in a chaotic attractor 
of the unperturbed system. By selecting an unstable periodic orbit, which has a desired transport 
property, and stabilizing it via the extended delay feedback method, we can control transport 
properties of the deterministic inertia ratchet system. Also, we show that the extended delay 
feedback method can be utilized for separation of particles in the deterministic inertia ratchet 
system as a particle's initial condition varies. 
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I. INTRODUCTION 

The ratchet effect, i.e., a directional motion of a par- 
ticle using unbiased fluctuations, has attracted much at- 
tention in recent years [J, Q ■ An early motivation in this 
field is to explain an underlying mechanism of molecular 
motors which transport molecules in the absence of ap- 
propriate potential and thermal gradients [|[ . Lately, the 
ratchet effect has been studied theoretically and exper- 
imentally in many different holds of science, e.g., asym- 
metric superconducting quantum interference devices Q , 
quantum Brownian motion [5|, Josephson-junction ar- 
rays @|, application for separation of particles Q, etc. 
It has been known that two conditions should be met to 
obtain the ratchet effect First, a system has to be 
in a non-equilibrium state by a correlated stochastic [|[ 
or a deterministic perturbation 0] . Second, the breaking 
of the spatial inversion symmetry is required. In doing 
so, an asymmetric periodic potential, named the "ratchet 
potential", is introduced. 

Recently, several works concerning the control of 
ratchet dynamics have been presented. The applying of 
a weak subharmonic driving in a deterministic inertia 
ratchet system was used to enlarge the parameter ranges 
where regular currents are observed [l(|, and the signal 
mixing of two driving forces was considered to control 
transport properties in a overdamped ratchet system [Til ] . 
Also, the effect of time-delayed feedback on the ratchet 
system has been studied. The anticipated synchroniza- 
tion was observed in delay coupled inertia ratchet sys- 
tems .12] and the stabilization of chaotic current to low- 
period orbits was presented, using time-delayed feedback 
methods, in the deterministic inertia ratchet system [l3l ]. 
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Starting with the work of Ott, Grebogi, and Yorke [T3 |. 
various methods for controlling chaotic dynamics have 
been developed [HI]. Particularly, Pyragas proposed a 
simple and efficient method, which utilizes a control sig- 
nal with a difference between the present state of the sys- 
tem and the previous state delayed by the period of an 
unstable periodic orbit (UPO) [la ]. This method is non- 
invasive in the sense that the control signal vanishes when 
the targeted UPO embedded in a chaotic attractor is sta- 
bilized. Some limitations on the Pyragas method have 
been reported [ijj and the modifications of the Pyra- 
gas method have been proposed to improve its efficiency 
[l8l [l|| . Socolar et al. presented a method utilizing in- 
formation from many previous states of the system, and 
this method is called the extended time- delay autosyn- 
chronization or the extended delay feedback (EDF) [181 ] . 
The stability and analytical properties of a delayed feed- 
back system have been investigated [2(| • 



In this paper, we aim to control transport properties 
of the deterministic inertia ratchet system. For this pur- 
pose, we have controlled a chaotic current of the system 
to a regular current by stabilizing an unstable periodic 
orbit which has a desired mean velocity, via the EDF 
method. Also, we have shown that the EDF method can 
be utilized for separation of particles in the determinis- 
tic inertia ratchet system as a particle's initial condition 
varies. The rest of the paper is organized as follows. In 
Sec. II, we have shown transport properties of the un- 
perturbed deterministic inertia ratchet system. In Sec. 
Ill, the system controlled via the EDF method has been 
presented and the linear stability analysis of a periodic 
orbit in the presence of the EDF has been considered. In 
Sec. IV, via the EDF method, we have shown achieve- 
ments of the desired transport properties of the system 
and a separation of particles has also been presented as 
varying particle's initial condition. The paper finishes 
with conclusions in Sec. V. 
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II. DETERMINISTIC RATCHET SYSTEM 

The deterministic inertia ratchet system is written as 
the following dimensionless equation, 



x + bx + V'(x) = acos(woi)- 



(1) 



Here, V'(x) denotes the derivative with respect to x and 
b is the friction coefficient, ujq and a are the frequency 
and amplitude of the driving force, respectively. Figure[T] 
shows an asymmetric periodic potential, i.e., the ratchet 
potential V(x) described by 

V{x) = C-( ^ (27r(a; - d)) + ^ sin (4ir(x - d)) ) /4tt 2 J, 

" (2) 
where d, S and C are introduced so that the ratchet po- 
tential has a minimum at x = with V(0) = 0. 

This system exhibits both regular and chaotic behav- 
iors, depending on parameters (a, b,LOo) [HI, H3, HH . In 
this paper, we vary only the parameter a, and set b = 0.1 
and wo = 0.67. General transport properties of the de- 
terministic inertia ratchet system are shown in Fig. O In 
Fig. HJa), we plot a bifurcation diagram of the strobo- 
scopic recording of a particle's velocity at t = kT, where 
k is a positive integer and T is the period of the driving 
force. The mean velocity of the system as a function of 
the parameter a is depicted in Fig. [DJb). As shown in 
Fig. [UJb), multiple current reversals occur as the am- 
plitude of the driving force is varied. It has also been 
observed that the current reversal is related to a bifur- 
cation from chaotic to regular regime [2l[ and that the 
type-I intermittency exists in this bifurcation [22| . 

When the system exhibits a regular behavior, the time 
required for a particle to move from one well of the po- 
tential to another is commensurable with the period of 
the driving force. Hence, the mean velocity of a regular 
current shows a locking phenomenon given as follows: 
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FIG. 1: Asymmetric periodic potential, i.e., the ratchet po- 
tential V(x) = C - (sm(27r(s - d)) + \ sin(47r(:r - d)))/4n 2 5 
with d = -0.19, 6 = sin(27r|d|) + sin(47r|d|), and C = 
-(sin(27rd) + 0.25 sin(47rd))/47r 2 <5. 
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FIG. 2: Bifurcation diagrams as a function of a at 6 = 0.1 and 
u>o ~ 0.67. In the region from a = 0.063 to 0.071, coexisting 
attractors are found; (a) the stroboscopic recording of particle 
velocity, (b) the mean velocity of current. 
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where L is the spatial period of the ratchet potential 
(L = 1, as shown in Fig. [T] then vi = j a ), T is the time 
period of the driving force, and — is an irreducible frac- 
tion (n, m G Z) [23(. vi is the fundamental locking veloc- 
ity corresponding to a particle's current which advances 
one well of the ratchet potential in a positive direction 
with the period of the driving force. As shown in Fig. 
[21 the system exhibits regular behaviors in some param- 
eter ranges; a period-1 orbit with (v) = (a = 0.06), 
a period-2 orbit with (v) — \vi (a = 0.074), a period-4 
orbit with (v) = —\vi (a = 0.081), and a period-2 orbit 
with (v) = — hvi (a = 0.092). When the system shows a 
chaotic behavior, the mean velocity of the chaotic current 
is almost zero averaged. 

It is worthy of note that there are various UPOs em- 
bedded in a chaotic attractor of the unperturbed system 
and that their mean velocities agree with Eq. ([3]). By 
stabilizing an UPO that has a desired mean velocity (i.e., 
written by specific n and m in Eq. ([3])), we can achieve 
a desired regular transport of the deterministic inertia 
ratchet system instead of the zero averaged chaotic cur- 
rent. In this system, the periodic orbit is defined by 
(x(t),v(t)) — (x(t + T),v(t + T)), where x(t) = x(t)(mod 
1) and r is the time period of the orbit. We are inter- 
ested in stabilizing some targeted UPOs among various 
UPOs that agree with Eq. ([JJ. In Fig. [3J we have shown 
several UPOs, which are desired to be stabilized, located 
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FIG. 3: (Color online) Unstable periodic orbits, period-1 
(black line), period-2 (blue line), period-3 (red line), and 
period-4 orbits (dotted black line) obtained from the unper- 
turbed system at a = 0.083. 



by the Newton method. For period-n orbits, we consider 
two UPOs that have the same period r = nT = 2n7r/wo 
with different mean velocities, where n = 1, 2, 3, 4: one is 
a positive current with the mean velocity (v) = -^vi and 
the other is a negative current with the mean velocity 
(v) = ——vi. Particularly, period-1 orbits include an os- 
cillating orbit confined in one well of the potential with 
the mean velocity (v) = 0. By selecting a specific UPO 
and stabilizing it via the EDF method, we can easily 
control transport properties of the deterministic inertia 
ratchet system. 



III. LINEAR STABILITY ANALYSIS OF 
PERIODIC ORBITS 

The deterministic inertia ratchet system controlled by 
the EDF method is described as 



x + bx + V'(x) = a cos(woi) + F, 



(4) 



where F is a control signal, i.e., the delayed feedback de- 
scribed by the particle's present velocity and the previous 
velocities delayed by multiples of the period of UPO. F 
is denoted by 



F = K({1 - R) R m ^x{t - mr) - x(t)) , (5) 



deviation from the periodic orbit £o(t) be 5£(t) — — 
£o(t). According to the Floquet theory, 5£(t) can be 
described as 



JV 



k=l 



(0) 



where Xk + iu>k is the Floquet exponent and ttfc(t) = 
Uk(t+r) is an eigenvector. is a constant and N is the 
dimension of the Poincare surface. For one such mode, 
one can obtain the following deviation relation (dropping 
the index k) . After the period r of the periodic orbit has 
passed, the deviation is described as 



6£(t + t) = exp ((X + iuj)(t + T))u(t + 

= exp((\ + iw)T)6$(t) = (A- 



(7) 



ifi)^(t), 



where A + ifl is the Floquet multiplier. When the de- 
lay terms are included, the phase space of the system 
becomes infinite dimensional and the system has an infi- 
nite number of Floquet multipliers. If the largest Floquet 
multiplier satisfies |Ai + < 1, i.e., the leading Flo- 
quet exponent Ai (Ai = In |Ai +iOi|/r) is less than zero, 
thereby the targeted UPO is stabilized. 
The time evolution of 5£(t) is given by 



8£ = 



1 

-V"(x) -b 



(8) 



° Q I ) K((l -R)J2 R^SZit mr) - *£(t)), 



where the matrix in the first term in Eq. 1(5)) is the Ja- 
cobian of the unperturbed system and the second term 
comes from the presence of the EDF. The delayed terms 
in Eq. (JS} can be eliminated and consequently the time 
evolution of the small deviation from the periodic orbit 
could be governed by 





-V"(x) 



k ] A t'" 



6£(jt) = A5Z{t), (9) 



where the ratio of geometric series has a constraint 
| A ^q| < 1 for convergence. For an elimination of the 
delay terms, we use the relation 



where A" is a strength of feedback, r is a delay time, which 
coincides with the period of the targeted UPO, and R 
(0 < R < 1) is a parameter that adjusts the distribution 
of each term's magnitude in the control signal. When 
R = 0, the EDF method is the same as the Pyragas 
method, i.e., F = K~(x(t — t) — x(t)) . 

Now, let us consider the linear stability analysis of a 
periodic orbit in the presence of the EDF. Let the small 



5£(t - nr) = (A + ifl)- n 6£(t), (n = 1, 2, 3, • ■ ■)• (10) 

Eq. ([9]) requires information of the targeted UPO. Hence, 
the Floquet multiplier is related to the eigenvalue prob- 
lem of the monodromy matrix $ r , which satisfies 



$ 4 = A$ t , $o = I- 



(11) 
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FIG. 4: The leading Floquet exponents of period-1 orbits; 
(a), (b), and (c) exhibit the leading Floquet exponents for the 
positive, the negative, and the confined currents as a function 
of K for the given R, respectively. 



The eigenvalue of <£> T defines the Floquet multiplier as 
follows: 
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FIG. 5: The leading Floquet exponents of period-2 orbits; (a) 
and (b) exhibit the leading Floquet exponents for the positive 
and the negative currents as a function of K for the given R, 
respectively. 
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FIG. 6: The leading Floquet exponent of the period-3 positive 
orbit as a function of K for the given R. 



det[$ T - (A + ifi)I] = 0. (12) 

The Floquet multipliers are obtained by numerically 
solving Eqs. ([5|).([TT]). and fT2"]l, The results of leading 
Floquet exponents are shown in Figs. [U [5j [H and [7] for 
period-1 (t = T = 2tt/ujo), period-2 (t = 2T), period- 
3 (t = 3T), and period-4 orbits (r = AT), respectively. 
Each figure shows the leading Floquet exponent Ai as 
a function of the strength of feedback K for different 
values of the control parameter R. The results tell us 



the stabilized region of control parameters (K, R, r) , in 
which the targeted UPO is stabilized (Ai < 0). 

As shown in Fig. [U the period-1 positive and the nega- 
tive currents cannot be stabilized by the Pyragas method 
(R = 0). The positive current is stabilized at R > 0.6, 
the negative current is stabilized at R > 0.4, and the 
confined current, which only exists for the case of the 
period-1 orbits, can be stabilized even at R — 0. Such a 
phenomenon is roughly explained by the different degrees 
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FIG. 7: The leading Floquet exponents of period-4 orbits; (a) 
and (b) exhibit the leading Floquet exponents for the positive 
and the negative currents as a function of K for the given R, 
respectively. 



of instability of UPOs in the unperturbed system (see Ai 
values at K = in Fig. |4j. With the larger degree of 
instability in the unperturbed system, the UPO can be 
stabilized with a larger R. The stabilized region of K for 
the given R increases as R increases and if R is not very 
large, then the minimum of Ai is deeper as R increases. 
The EDF method (0 < R < 1) is more effective than the 
Pyragas method (R = 0). 

Another interesting property of the EDF method is 
a limitation on the minimum of Ai for given control 
parameters R and r. As shown in Fig. [4J the lead- 
ing Floquet exponent is always larger than A*, where 
A* = ln(i?)/r. A* is given by the constraint on the ratio 
of geometric series. In Eqs. (|5j)- (fTU)) , the system gains 
the limit of geometric summation of the delayed feedback 
terms, which do not diverge, when the Floquet exponent 
is infinitesimally larger than A*. A* of the period-I or- 
bits (r = T = 2tt/ujo) is numerated as A* ~ —0.1716 
(R = 0.2), -0.0977 (R = 0.4), -0.0544 (R = 0.6), and 
-0.0237 (R = 0.8). If the result of the Floquet exponent 
obtained by solving Eqs. (O, (fTTj) , and (fT2")) is far from 
A* for given R and r, then the shape of Ai(AT) makes a 
narrow valley. If the result is very close to A*, then the 
shape of \\{K) makes a valley with a flat bottom. In 
this case, the minimun of Ai is infinitesimally larger than 
A* and Ai is very slowly increasing as K increases. 




FIG. 8: The loci of the largest Floquet multiplier of period- 
2 negative orbit as K varies from to oo; (a) loci at R = 
0.2 and (b) loci at R = 0.6. The outer (inner) circle has a 
radius of 1 (R). The open circle and the crosses denote the 
location of the largest Floquet multiplier at K = and the 
discontinuity, respectively. For K — » oo, the largest Floquet 
multiplier approaches (Ai,fii) = (1,0) in both cases of (a) 
and (b). 



The leading Floquet exponents for period-2, period- 
3, and period-4 orbits are shown in Figs. OH and 
respectively. The results are qualitatively equivalent to 
those of period- 1 orbits shown in Fig. 2] The period- 
2 positive current is stabilized at R > 0.4, and the 
period-2 negative current is stabilized at R > 0.2. A* 
of the period-2 orbits (r = 2T = Att/ujo) is numerated 
as A* ~ -0.0858 (R = 0.2), -0.0488 (R = 0.4), -0.0272 
(R = 0.6), and -0.0118 (R = 0.8). The period-3 nega- 
tive current is stabilized at R > 0.8. A* of the period-3 
orbit (r = 3T = 6tt/ujq) is numerated as A* ~ —0.0572 
(R = 0.2), -0.0325 (R = 0.4), -0.0181 (R = 0.6), and 
—0.0079 (R = 0.8). The period-3 positive current cannot 
be stabilized by the EDF method because the largest Flo- 
quet multiplier is a positive real value when the system 
is unperturbed. It has been known that the delayed feed- 
back method, including the EDF method, can stabilize 
only a certain class of periodic orbits with a finite torsion 
(171 ]. The period-4 positive current cannot be stabilized 
at R < 0.8 because it has a greater degree of instabil- 
ity (Ai ~ -862, Qi = 0) when K = 0. The period-4 
negative current can be stabilized even at R = 0. A* of 
the period-4 orbits (r = AT = 8it/ujo) is numerated as 
A* ~ -0.0429 (R = 0.2), -0.0244 (R = 0.4), -0.0136 
(R = 0.6), and -0.0059 (R = 0.8). 

All periodic orbits in Figs. and[7]have common 

properties. When K = 0, their largest Floquet multipli- 
ers are real and negative so that the corresponding lead- 
ing Floquet exponents satisfy u)\ = tt/t. It means that all 
UPOs flip their neighborhood within the period r in the 
unperturbed system. The variation of each orbit's largest 
Floquet multiplier depending on the feedback strength K 
has a common aspect as shown in Fig. [5J We plot the 
loci of the largest Floquet multiplier of period-2 nega- 
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FIG. 9: Stabilized period- 1 orbits (insets) and the dynam- 
ics of the control signal F at K = 0.67, R = 0.8, and 
t — T = 2n/uio', (a) the positive current from the initial con- 
dition (xo,vo) = (—0.35,0.2), (b) the negative current from 
(0.2,0.0), and (c) the confined current from (0.0,0.0). 



tive orbit at R = 0.2 (Fig. Ufa)), and at R = 0.6 (Fig. 
HJb)). Let us see Fig. Ufa). As K increases, the largest 
Floquct multiplier moves toward the zero point in re- 
maining a negative real value and crosses the unit circle 
so that the period-2 negative orbit is stabilized (Ai < 0). 
With the further increase of K, the largest Floquet mul- 
tiplier becomes a complex value. It is precisely at this 
point that the leading Floquet exponent Ai has a mini- 
mum value (see Fig. [5jb) for the case of R — 0.2). Then 
the pair of complex conjugates move symmetrically and 
walk out of the unit circle so that the periodic orbit is 
destabilized. For K — > oo, the largest Floquet multiplier 
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FIG. 10: Stabilized period-2 orbits (insets) and the dynam- 
ics of the control signal F at K = 0.35, R = 0.8, and 
r = 2T = Ak/ojo; (a) the positive current from the initial 
condition (xo,vo) = (0.0,0.0), and (b) the negative current 
from (0.0,0.1). 



approaches (Ai,fii) = (1,0) very slowly. In the same 
manner with Fig. [HJa), the largest Floquct multiplier in 
Fig. [Spa) (at R — 0.6) moves toward the zero point as 
K increases, until it meets the inner circle of radius R. 
Then, it jumps to the opposite side of the inner circle and 
becomes a complex value. In some finite interval of K, 
the pair of complex conjugates follow a curved line, which 
is very close to the inner circle. The loci of the largest 
Floquet multiplier in this interval explain the presence 
of a flat bottom in Fig. [5fb) for the case of R = 0.6. 
(When the Floquet multiplier locates on the inner cir- 
cle of radius i?, the corresponding Floquet exponent is 
A = In |A+iO|/r = \n(R)/r = A*.) Finally, the pair walk 
out of the unit circle and approach (Ai,fii) = (1,0). 



IV. CONTROL OF TRANSPORT PROPERTIES 

With the results of linear stability analysis of periodic 
orbits, we can obtain a desired transport property of the 
system (i.e., a regular current with a desired mean veloc- 
ity) by choosing the control parameters (K, R, r) where 
the corresponding UPO is stabilized. For simple and ef- 
ficient application of the EDF method, it is hoped that 
each UPO has its own stabilized region of control param- 
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FIG. 11: (Color online) (a) Basins of period-1 orbits; the ini- 
tial points marked by □ (blue), A (red), and • (black) are the 
basins of the positive, the negative, and the confined currents, 
respectively, (b) Basins of period-2 orbits; the initial points 
marked by ■ (black) and A (red) are the basins of the positive 
and the negative currents, respectively. 



eters, in which the other UPOs still remain in an unstable 
state. Some of the UPOs have their own stabilized re- 
gions of control parameters. These orbits are the period- 
1 confined, the period-2 negative, the period-3 positive, 
and the period-4 negative currents. For the cases of the 
period-1 positive, the period-1 negative, and the period-2 
positive currents, the stabilization of each periodic orbit 
is rather complex because they do not have their own 
stabilized regions of control parameters. The stabilized 
region of the period-1 positive and the negative currents 
always overlaps with that of the period-1 confined cur- 
rent, and the stabilized region of the period-2 positive 
current overlaps with that of the period-2 negative cur- 
rent. 

Now, we are interested in the multistable phenomenon 
that more than one UPOs are stabilized at the same 
control parameters (K,R,t). At control parameters, 
K = 0.67, R = 0.8, and r = T = 27r/w , all of period- 
1 orbits are stabilized. In Fig. 2J each of the leading 
Floquet exponents of the period-1 orbits at theses pa- 
rameters is marked by a black circle and all of them are 
less than zero. Also, all of period-2 orbits are stabilized 
at K = 0.35, R = 0.8, and r = 2T = Att/loq. Each of the 
leading Floquet exponents of the period-2 orbits at these 



parameters is marked by a black circle in Fig. [5l In the 
unperturbed system {K = 0), all of the initial conditions 
evolve into a chaotic current in the same manner. How- 
ever, the system controlled by the EDF method shows 
different currents for different initial conditions (xq,Vq), 
at control parameters where the multistable phenomenon 
is observed. 

Before considering the numerical integration for ob- 
taining dynamics of the system controlled by the EDF 
method, we rewrite the control signal F(t) given in Eq. 
([3]) into a more convenient form 

F(t) = K((l-R)S(t-T)-x(t)), (13) 
S(t) = x(t) + RS{t-T), 

where S(t) = X)m=o R m x{t — mT) for an equivalent equa- 
tion with Eq. ([5]) (see Ref. [24|). In the following numer- 
ical integrations, we set S(t ) — for t in the interval 
[— r, 0] and initialize S(t ) = x(t )/(l — R) for t in the 
interval [0, r]. Then, the system is not perturbed (F = 0) 
for t in the interval [0, r] and perturbed by the control 
signal from t = r. In Fig. [9l we have plotted three sta- 
bilized period-1 orbits that evolved from different initial 
conditions and the dynamics of the control signal F at 
the same control parameters, K = 0.67, R — 0.8, and 
t = T = 2tt/luq. For each initial condition, we have in- 
tegrated Eqs. (g]) and (fT3)l . In Fig. [TU1 we have shown 
two stabilized period-2 orbits that evolved from different 
initial conditions and the dynamics of the control signal 
F at K = 0.35, R = 0.8, and r = 2T = 4tt/w . The mul- 
tistable phenomenon shows that the EDF method can 
be utilized for separation of particles in the determinis- 
tic inertia ratchet system. Via the EDF method, we can 
separate particles in the deterministic inertia ratchet sys- 
tem as their initial conditions vary. In Fig. 111! we have 
investigated the basins of period-1 and period-2 orbits. 
We have integrated Eqs. ([4]) and (fT3|) with the initial 
condition (xq,Vq), where xq is distributed in one well 
of the potential, x £ (—0.38,0.62), and v Q is confined 
to the ranges of velocities in the unperturbed system, 
v e (—0.3,0.3). As shown in Fig. [TlTa). the basins 
of the period-1 positive, the negative, and the confined 
currents are marked by □ (blue), A (red), and • (black), 
respectively. The basins of the period-2 positive and the 
negative currents are marked by • (black) and A (red), 
respectively, in Fig. [TTT b). 

V. CONCLUSIONS 

We have studied the control of transport properties in 
the deterministic inertia ratchet system via the extended 
delay feedback method. We have controlled a chaotic 
current of the unperturbed system to a regular current, 
which has a desired mean velocity. To obtain the control 
parameters in which the corresponding unstable periodic 
orbit is stabilized, we solve the leading Floquet exponent 
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in the presence of the extended delay feedback. With 
the results of leading Floquet exponents as a function of 
control parameters, we have obtained a desired regular 
transport property of the system. Also, we have observed 
the multistable phenomenon that more than one unstable 
periodic orbits are stabilized at the same control param- 
eters and we have shown that the extended delay feed- 
back method can be utilized for separation of particles 
as a particle's initial condition varies. 
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